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We present a mechanism for generating primordial magnetic fields with large correlation lengths 
on the order of 100 kpc today. The mechanism is based on recently conjectured QCD domain 
walls or similar CP violating domain walls with QCD scale structure. Such domain walls align the 
electric and magnetic dipole moments of the nucleons shortly after the QCD phase transition, pro- 
ducing electromagnetic fields correlated along the domain walls. Through the Kibble mechanism, 
the domain walls attain Hubble-scale correlations which they transfer to the aligned electromagnetic 
fields. Due to the CP violation, the Hubble-scale walls produce helical (non-zero Chern-Simons) 
magnetohydrodynamic turbulence which then undergoes an inverse cascade, allowing the correla- 
tion lengths to grow to 100 kpc today. We present an estimate the magnitude of the generated 
electromagnetic fields in terms of the QCD parameters. We also discuss some other unexplained 
astrophysical phenomena which may be related to this mechanism. In particular, we discuss the 
relation between primordial magnetic fields and the Greisen-Zatsepin-Kuzmin (GZK) cutoff viola- 
tions. We also outline how, by creating inhomogeneities in the nucleon density, QCD domain walls 
may lead to inhomogeneous big bang nucleosynthesis (IBBN) explaining the Q.b excess recently 
measured through cosmic microwave background (CMB) distortions. 
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I. INTRODUCTION 

The purpose of this paper is to present a careful ar- 
gument of the mechanism outlined in our letter to 
generate large scale cosmic magnetic fields. It is an ex- 
tended and updated version of the letter |l| and the con- 
ference proceedings |^ and should be considered to be 
the definitive version of these papers. 

Many different observations suggest that there exist 
substantial (microgauss) magnetic fields in the universe 
today H, however, there has yet to emerge a the- 
ory which adequately explains the origins of these fields. 
Most of the data on large scale astrophysical magnetic 
fields comes from the observation of synchrotron radia- 
tion emitted in galaxies. This radiation is plane polar- 
ized, and as it passes through magnetic fields, the plane 
of polarization rotates due to the Faraday effect: an ef- 
fect which depends on the frequency of the radiation and 
the strength and orientation of the magnetic fields. By 
comparing several sources or radiation with different fre- 
quencies, one can extrapolate to determine the original 
plane of polarization and then estimate the magnetic field 
strengths. 

What is striking, is not just the existence of magnetic 
fields, but that they appear to be microgauss fields which 
have correlations as large as 500 kpc in clusters. To put 
this into perspective, the luminous core of galaxies have 
typical scales of up to 10 kpc while it is estimated that 
the galactic dark matter halos extend to 50 kpc. Thus, it 
seems that galactic dynamo mechanisms cannot produce 
these large correlation lengths. 

The current models for producing these fields involve 
two main processes: 1) dynamical amplification and/or 
generation of magnetic fields by galactic processes (galac- 
tic mechanisms) and 2) primordial mechanisms which 



take place prior to gravitational structure formation. The 
galactic mechanisms are primarily based on gravitational 
dynamos, although there are suggestions that supernovae 
or other stellar phenomena may play a role. While it is 
likely that galactic dynamos amplify fields, it seems dif- 
ficult to account for the large scale correlations of the 
magnetic fields when only galactic mechanisms are con- 
sidered. It is also not certain that galactic mechanisms 
can generate magnetic fields: instead they serve only as 
an amplifier requiring seed fields to be present for the 
dynamo to work. 

The inadequacies of the galactic mechanisms have lead 
to many proposals that the magnetic fields may have a 
primordial origin. In this case, some process in the early 
universe (typically at a cosmic phase transition or dur- 
ing infiation) is thought to generate magnetic turbulence. 
This turbulence then sustains itself as the universe ex- 
pands and what we observe are the remnants of this tur- 
bulence. Most primordial sources, however, also produce 
fields which end up with very small correlations today or 
which are very weak. 

Most likely, a complete picture of the history of astro- 
physical magnetic fields requires some primordial inputs 
as well as amplification through gravitational dynamics. 
In this paper, however, we discuss a primordial mecha- 
nism which seems to naturally produce fields of 100 kpc 
correlations today. In combination with dynamic amplifi- 
cation mechanisms, we hope that this mechanism might 
provide a solid foundation for the theory of large-scale 
astrophysical magnetic fields. 

At this point we would like to refer the reader to the 
several reviews and sample papers in this field. The pri- 
mary discussions of observations, which contain reviews 
of the theory, are presented in |^. Good current re- 
views are given in ||, ^. Many different types of 
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primordial mechanisms are discussed, for example: In- 
flationary mechanisms H, cosmic strings jgl ITol, c harge 
asymmetries |2j , and phase transitions |13|, ij] . The 
evolution of primordial magnetic fields is discussed in 
[0, |l^, |l6|. In particular, the inverse cascades discussed 
here will be important for our mechanism. 



II. OVERVIEW OF THE MECHANISM 

The mechanism that we propose has the following core 
components: 

1. Sometime near the QCD phase transition, Tqcd ~ 
1 GeV, domain walls form which can interact with 
QCD scale physics. 

2. These domain walls rapidly coalesce until there re- 
mains, on average, one domain wall per Hubble vol- 
ume with Hubble scale correlations. 

3. Baryons interact with the domain walls and align 
their spins along the walls. The Hubble scale cor- 
relation of the domain walls thus induces a Hubble 
scale correlation in the spin density. 

4. The anomalous magnetic and induced electric 
dipole moments of the baryons generate helical 
electromagnetic fields also correlated on the Hubble 
scale. 

5. The domain walls move rapidly and vibrate, effec- 
tively filling the Hubble volume with helical mag- 
netic turbulence with a Hubble scale correlation. 

6. The domain walls decay and the electric fields are 
screened leaving magnetic turbulence with Hubble 
scale correlations. 

7. As the universe expands, an "inverse cascade" 
mechanism transfers energy from small to large 
scale modes, effectively increasing the resulting cor- 
relation lengths but diluting the field strengths. 

8. Galactic dynamos amplify the fields in galaxies, but 
the fields should also persist in the extra-galactic 
media. 

The idea that domain walls might generate magnetic 
fields is not original. For example, it was suggested 
that standard axion domain walls could be ferromagnetic 
in jTrt , however mechanism discussed their seems to be 
flawed: The scale of the standard axion walls is of the 
order which is at least some twelve orders of magni- 
tude larger than the QCD scale AgJ^p {rua is the axion 
mass) . It is hard to see how these walls can efficiently af- 
fect QCD physics at the temperatures that were present 



in the early universe where thermal fluctuations will de- 
stroy all coherence.^ 

Another problem with proposals including standard 
axion domain walls is that these walls must decay to pre- 
vent cosmological problems . There are still questions 
about how the standard so-called N ^ 1 axion domain 
walls can decay. 

In this paper, we outline the properties that domain 
walls must have to generate sufficient magnetic seed 
fields. The exact source could be one of several types 
of walls, including modified axion domain walls, or en- 
tirely different types of walls. To be concrete, we present 
our model in terms of a recently conjectured quasi-stable 
QCD domain wall [po[ which may exist, either indepen- 
dently of axion physics, or along with a dynamical axion 
adding additional QCD scale structure to the standard 
axion domain walls. These domain walls are character- 
ized by a transition in the singlet rj' field which has a 
size and energy scale set by Aqcd- Hence, QCD domain 
walls can directly couple to QCD physics. In addition, 
the singlet field transition at the center of the wall in- 
duces an effective non-zero CP violating 9 background 
which in turn will induce an electric dipole moment and 
alter the magnetic dipole moment in the fermions |2l| so 
that both the electric and magnetic dipole moments of 
all the particles are on the same order. In the presence of 
these anomalous dipole moments, the cancellations dis- 
cussed in p^ , |23[ | in the domain wall background are no 
longer a problem. 

Another crucial aspect of our mechanism is the "in- 
verse cascade" which governs the evolution of the mag- 
netic fields after they are formed. This mechanism was 
suggested by Cornwall Q, discussed by Son and 
confirmed by Field and Carroll |2^. It is based on the 
idea that magnetic helicity (Abelian Chern-Simons num- 
ber) H — / A • Bd'^x is approximately conserved in the 
universe where temperatures are higher than Tq « 100 
eV. This conservation of helicity causes energy to cas- 
cade up the turbulent modes increasing the energy in 
large scale modes and increasing the effective correla- 
tion length of the turbulence. The importance of helicity 
was originally demonstrated by Pouquet and collabora- 
tors Without this helical inverse cascade, there is 
no known way to generate large correlations fields today 
from sub-Hubble scale fields in the early phase transition 
and one must consider super-Hubble scale correlations 
resulting form inflationary schemes. It turns out, how- 
ever, that Hubble scale correlations at the QCD phase 
transition (the last major phase transition) provide cor- 
relations on the order of 100 kpc today: thus it is natu- 
ral to consider QCD physics as the source of primordial 



^ For some other discussions about the magnetic properties of the 
domain walls, see ||l8|| and reference therein. We should note, 
however, that in all these discussions, the most difficult problem 
of generating large scale correlations has not been addressed. 
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fields (earlier physics can only produce even smaller cor- 
relations). 

Perhaps the least understood aspect of this mechanism 
concerns the dynamics of the domain walls and the in- 
teractions of the domain walls, nucleons and electric and 
magnetic fields. As we shall show, all of these compo- 
nents interact on the same scale of Aqcd and hence there 
are complication back-reactions and nonlinear dynamics. 
We presently do not have the tools to fully analyze these 
features, but we present here in detail quantitative esti- 
mates and calculations which we believe are good esti- 
mates of the scale of the effects. 

The result is a mechanism which naturally produces 
magnetic fields today with / ~ 100 kpc correlations and 
with strengths of 



lO-^G 

^Aqcd ' 



I ~ 100 kpc 



(1) 



where the parameter 1 < ^Aqcd 10^^ depends on 
the dynamics of of the domain walls as discussed in sec- 
tion (IV). If the correlation ^ turns out to be small, then 
this mechanism might generate detectable extra-galactic 
fields, otherwise we still require a galactic dynamo to 
amplify the fields. In any case, however, these seeds 
still maintain the large scale correlations of the observed 
fields, and it seems that even if ^ is large, the resulting 
fields may be strong enough to seed the galactic dynamos 

We shall begin by discussing the inverse cascade mech- 
anism in section (III) and then give estimates of the field 
strengths in an idealized case of static, flat walls. Finally, 
we shall discuss the dynamics of the domain walls and 
describe the whole process, justifying the mechanism. 



III. EVOLUTION OF MAGNETIC FIELDS 

Given a stable magnetic field configuration in the uni- 
verse, one might naively expect the size of the correla- 
tions of the field to expand with space as governed by 
the universe's scale parameter / cx R{T) and the field 
strength to be correspondingly diluted B cx R{T)^'^. It 
was discovered by Pouquet and collaborators p7| , how- 
ever, that if the magnetic fields have a non-zero helic- 
ity (Abelian Chern-Simons number) H — / A • Bd'^x, 
then the fields will scale differently. Cornwall |24| sug- 
gested that helical fields might undergo an inverse cas- 
cade. The magnetohydrodynamic (MHD) equations were 
subsequently studied by Son who derived the scaling 
relations (^) presented below. These have subsequently 
been confirmed by Field and Carroll . 

The reason for the scaling is that the magnetic helicity 
H is an approximately conserved quantity in the early 
universe. It is also known that the small scale turbulent 
modes decay more rapidly than the higher scale modes. 
In order to conserve the helicity as the small scale modes 
decay, the helicity must be transfered to larger modes 



and with this transfer of helicity is a transfer of energy. 
This is the source of the inverse cascade. 

To understand the origins of the conservation of helic- 
ity, note that it is a topological quantity that describes 
the Gaussian linking number of the vector potential lines 
of flux (see for example [^). In a perfectly conducting 
medium, these lines of flux cannot cross, and hence there 
is no way to unlink the flux lines: helicity is perfectly con- 
served. Even when the conductivity is flnite, the helicity 
is also well conserved. 

There is a direct analogy with fluid mechanics. The 
equivalent there is vorticity C = / v • (V x v)d'^a; which is 
the Gaussian linking number of the fluid flow lines. If the 
fluid has no viscosity, then ( is perfectly conserved. Even 
viscous fluids, however, approximately conserve vorticity. 
This is why, for example, smoke rings and tornados are 
so stable. 

We shall not derive the scaling relationships here, in- 
stead we refer the reader to |2^, |2^, however, we sum- 
marize here the results: 
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(2b) 



These relate the initial field strength -Brms(Ti) with ini- 
tial correlation l{Ti) to the present fields today (Tnow ~ 
2 X 10^"* cV) -Brms(Tnow) with Correlation ^(Tnow)- Dur- 
ing the period when the universe supports turbulence 
(as indicated by a large Reynolds number Re), the in- 
verse cascade mechanism functions and have the scalings 
B cx T''/'^ and I oc T~^/^ as indicated by the second fac- 
tors in (|^). In the early universe. Re is very large and the 
turbulence is well supported. As the universe cools, even- 
tually, for temperatures below some Tq, the turbulence is 
no longer well supported. Exactly what the effective tem- 
perature To the turbulence ceases is not clear: Son points 
out that at Tq ~ 100 cV, the Reynolds number drops to 
unity and thus turbulence is not well supported because 
of the viscosity of the plasma |2^. We take this as a 
conservative estimate. Field and Carroll argue that the 
turbulence is force-free and so unaffected by the viscosity. 
Thus they take Tq ~ 1 eV at the epoch when the matter 
and radiation energy densities are in equilibrium and ar- 
gue that the cascade may even continue into the matter 
dominated phase of the universe. If this is true, then it 
might be possible to increase the correlation lengths of 
the fields by one or two orders of magnitude from the 
conservative estimate (0). In any case, for temperatures 
lower than Tq, the turbulence and inverse cascade are not 
supported and so we assume that the fields are "frozen 
in" and experience only the nai've scaling I cx and 
B ocT^ indicated by the first factors in (0). 

As pointed out by Son the only way to gener- 
ate turbulence is either through a phase transition Tj 
or through gravitational instabilities. Thus, until grav- 
itational dynamos are active, the scalings (0) should be 
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valid. In any case, galactic dynamos will amplify the 
fields, but will not affect the correlation length, so in 
particular, (|^) should be a good estimate, regardless of 
galactic dynamics (the uncertainty come during the tran- 
sition period T ^ when the scaling laws change. The 
estimate (|^) is generally considered conservative in this 
sense as the cascade likely continues for some time past 
To @]. 

Now we consider the source of the magnetic turbulence 
at a phase transition. As we shall show, our mechanism 
generates Hubble size correlations U at a phase transition 
Tj. In the radiation dominated epoch, the Hubble size 
scales as T~^. Combining this with ([2b|), we see that 



^now ^ T- 



-1/3. 



thus, the earlier the phase transition, the 



smaller the resultant correlations. 

The last phase transition is the QCD transition, 
Ti = Tqcd ~ 0.2 GeV with Hubble size ^(rQCD) ~ 
30 km. With our estimates (|4^) of the initial mag- 



netic field strength Br^siTi, 



eA 



<CD 



/(^Aqcd) 



(10^''G)/(^Aqcd) we use Equations (||) to arrive at the 
estimate (^. The meaning of the correlati on le ngth f 
will be discussed in detail later in section (VII). The 
most important result here is that, as long as one has a 
mechanism to generate Hubble scale correlations and a 
maximally helical magnetic field at the QCD phase tran- 
sition, magnetic turbulence of lOO kpc correlations is nat- 
urally produced. The questions: 'How can helical mag- 
netic fields with Hubble-scale correlations be produced at 
the QCD phase transition' and, 'Are these fields strong 
enough to account for the observed microgauss fields?' 
will be addressed in the rest of this paper. The estimate 
(|I|) suggests, however, that even in the worst case of al- 
most maximal suppression ^Aqcd ~ 10"'^^, an efficient 
galaxy dynamo may be able to amplify the fields to the 
microgauss level. In the best case, the mechanism would 
produce measurable extra-galactic fields. 

In either case, the important result is the generation 
of the lOO kpc correlations: if observations show that 
the fields have much larger correlations, then the pro- 
posed mechanism can only be salvaged if a more efficient 
"inverse cascade" mechanism is shown to work between 
Tqcb and now. Having said this, one might consider 
the electroweak or earlier phase transitions. As we men- 
tioned, the earlier the phase transition, the smaller the 

— 1/3 

resulting correlations Inow . For the electroweak 

transition, the scaling ( pb| ) suggests that Hubble scale 
helical fields could generate 100 pc correlations today. 
Thus it might be possible that electroweak phenomena 
could act as the primordial source, but this presupposes 
a mechanism for generating fields with Hubble scale cor- 
relations. Such a mechanism does not appear to be pos- 
sible in the Standard Model. Instead, the fields produced 
are correlated at the scale which can produce only 
~ 1 km correlations today which are of little interest. 

Thus, the previous analysis seems to suggest that, in 
order to obtain magnetic fields with 100 kpc correlation 
lengths, helical fields must be generated with Hubble 
scale correlations near or slightly after the QCD phase 



transition Tqcd ■ The same conclusion regarding the rel- 
evance of the QCD scale for this problem was also reached 
by Son, Field and Carroll ||2|, |2^. Without further ado, 
we now present our picture of the mechanism and justify 
the the estimate (□) of the magnetic field strength. 



IV. GENERATION OF MAGNETIC FIELDS BY 
DOMAIN WALLS 



The key players in this mechanism are domain walls 
which form shortly after the QCD phase transition. De- 
tails of the walls were presented in |Q which wi ll be 
summarized in section (^). In sections (VI) and (VH) 
we shall show that these walls tend to align nuclear mag- 
netic and electric dipole moments along the plane of the 
wall. An important feature of the walls is that across 
the wall there is maximal strong CP violation due to 
an induced nonzero 6. Because of this, the electric and 
magnetic dipole moments of the nucleons are of the same 
order. Thus, both neutrons and protons will have non- 
zero electric and magnetic dipole moments and play a 
role in generating the electromagnetic fields. 

Because of the correlation between the electric and 
magnetic fields along the domain wall, the generated 
fields have an induced helicity as we shall examine in 
section (VH). This helicity has the same sign along the 



entire domain wall and we expect that the domain wall 
will fill the entire Hubble volume, thus the helicity will 
be correlated on the Hubble scale. 

Finally, the domain walls will decay as discussed in 
pot so that the universe is not dominated by domain 
walls today. By this point, however, the helical magnetic 
turbulence has been generated. 



A. Hubble Size Correlations 

The reason that we feel domain walls hold the key to 
explaining primordial magnetic seeds is that in a short 
time they can generate Hubble scale correlations. The 
initial fields must have a Hubble scale correlation or else 
there is no known way — even with the inverse cascade — 
to generate the huge correlations today. Let us briefiy 
summarize what we expect to be the behaviour of domain 
walls at the QCD phase transition. For a nice description 
of general domain wall dynamics see from which most 
of these results were derived. 

1. Prior to the phase transition Ti — Aqcd, the fields 
are in random fiuctuations on the scale Ti and do- 
main walls are not present. 

2. After the phase transition, however, the fields set- 
tle into their vacuum states. Domains are formed 
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where the fields are setthng into different^ vacuum 
states. These domains are separated by domain 
walls and have a scale set by Aqcd- 

3. Numerical studies suggest that these small-scale 
domain walls rapidly merge increasing the corre- 
lation length of the walls. This coarsening occurs 
simultaneously throughout space and the correla- 
tion length of the domain walls can increase faster 
than the speed of light. 

4. The coarsening stops once the domain walls attain 
a Hubble scale. On average, one ends up with one 
domain wall per Hubble volume, but which curls 
and moves through space, essentially filling the vol- 
ume. 



It is these Hubble sized domain walls that can gener- 
ate magnetic turbulence with Hubble size correlations. 
As we shall see below, there are two types of domain 
walls corresponding to opposite field transitions. One we 
call a "soliton" , and the other we call an "anti-soliton" . 
Together a soliton and an anti-soliton can annihilate, 
but the coarsening essentially separates regions of soli- 
tons from anti-soliton regions by a distance of the Hubble 
scale so that they do not annihilate. In section (VILA) 
we shall show that the solitons and anti-solitons are as- 
sociated with helicity of the opposite sign. Thus, the 
domain walls effectively separate the helicity generating 
a Hubble scale correlation length in the fields and in the 
helicity.'^ 



In the case of QCD domain walls, the vacuum states are actually 
the same but the field configuration going from one domain to the 
next undergoes a classically stable transition. This behaviour is 
qualitatively similar to the sine-Gordon model Csg = (9^0)^ — 
cos <f> where <f> is interpreted as a phase so that the vacuum states 
<f> = 2-Kn are actually identical. 

There is some question about what conditions must be like at the 
phase transition in order for domain walls to form. QCD lattice 
simulations suggest that at low densities (such as those present 
in the early universe), the transition between the quark- gluon 
plasma and the normal hadronic phase is a smooth crossover 
and that the critical point sits at some finite density (see the 
recent review PCj]). If the rate at which the universe cools is suf- 
ficiently slow, tnen it is possible that no domain walls will form. 
In the preface to the paperback edition [ p^ , the authors discuss 
this scenario as the Kibble-Zurek picture: to estimate the size 
of the correlations produced, one must consider the relaxation 
timescale t{T) that it takes to establish correlations on the scale 
^(r). The freezout temperature Tj is determined by the condi- 
tion T ~ ij5 = |T — Tc|/|T|, i.e. when the relaxation time is 
on the same order as the dynamical timescale of temperature 
variations. Since the true critical point is at a somewhat higher 
densities than the universe, tjj my be bounded from below and 
if the cooling is sufficiently slow, it is possible that r <C to and 
domain walls will not form. To estimate these effects requires a 
better understanding of the dynamics of the domain walls and 
of the phase transition than we presently have. We assume that 
the dynamics are such that domain walls do form and coalesce 
forming Hubble-scale correlations as described in |19|. 



B. Essential Domain Wall Properties 

Thus, we can formulate a set of properties that must 
be satisfied by domain walls if they are to be considered 
as sources for the primordial magnetic fields described in 
this paper: 

1. The walls must attain Hubble scale correlations 
near the QCD phase transition to generate the ob- 
served correlations. 

2. The walls must have structure on the scale of Aqcd 
in order to interact effectively with nucleons. 



3. 



There must be some way to avoid the cancellations 
discussed in 



23 so that the walls are ferromag- 



netic rather than diamagnetic. 

4. The walls must somehow induce a definite helicity 
throughout the Hubble volume. 

5. The domain walls must be unstable or have other 
features so that the problems of domain wall dom- 
ination in the universe are avoided, but they must 
be stable enough that they can generate the appro- 
priate fields. They must also leave nucleosynthesis 
production ratios relatively unaffected. 

It seems to be a rather general property of cosmological 
domain wall networks that they rapidly coarsen through 
the Kibble mechanism until the walls have a Hubble-scale 
correlation length Thus, criterion |l| should be easily 
satisfied by almost all types of domain walls. Criterion 
PI ru les out the standard axion domain walls discussed in 
|l7| , |3l] , however, there may be features of these walls 
that have QCD scale which were previously neglected. 
In particular, Shifman and Gabadadze discuss a gluonic 
core sandwiched at the center of axion domain walls ]3^ 
which they call D-walls. This structure has a QCD scale 
and may be able to align nuclear matter. 

In our paper pc| | we discuss another possibility: that 
the rj' field may provide domain wall structure with QCD 
scale. This structure, which we shall refer to as a QCD 
domain wall, can exist within the standard axion domain 
wall providing the required QCD scale structure, but can 
also exist, even if there is no axion (unlike the D-walls of 
p3| which require an axion field). 

A further property of QCD domain walls and axion do- 
main walls is that, at their center, there is strong CP vio- 
lation. This CP violation has several effects as described 



in section ( [ V C ) . In particular, the CP violation induces 
an anomalous electric dipole moments in nucleons. Thus, 
CP violating domain walls can satisfy criterion ^. 

Helicity is also associated with CP violation as it is a 
CP odd quantity. As discussed above, the soliton and 
anti-soliton domain walls solutions have opposite CP. 
Thus, each is associated with opposite helicity. Typi- 
cally, the solitons and anti-solitons separate spatially so 
that Hubble-sized regions are filled with one type or an- 
other. The helicity is generated through the correlation 
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of both electric and magnetic fields along the walls. Thus, 
the Hubble scale spatial separation of soliton and anti- 
soliton domain walls also separates the helicity and thus 
generates helical turbulence with Hubble scale correla- 
tions satisfying criterion ^. 

Another major problem with axion domain walls is 
that most varieties appear to be absolutely stable. The 
N = 1 axion model discussed recently by Chang, Hag- 
mann and Sikivie has a decay mode that satisfies the 
criterion p. However, this is not a phenomenologically ac- 
ceptable model and, to date, the other axion models are 
plagued by this problem. If axion domain walls are to be 
considered, then a satisfactory solution to this problem 
must be found. 

The QCD domain walls |^ without an axion, exhibits 
a similar decay mode to the the A'^ = 1 axion model. The 
scales, however, are set by Aqcd rather than ma and so 
criterion |^ must be address from the point of view: Do 
the walls live long enough to generate the turbulence? 
As addressed in Q, the answer may be yes. The only 
decay mode is through a nucleation process suppressed 
by quantum mechanical tunneling. Consequently, these 
walls may have a macroscopic lifetime long enough to 
generate the fields. In any case, however, they decay 
fast enough to avoid affecting nucleosynthesis and other 
cosmological effects. 

Thus, there may be several types of domain walls that 
could act as sources for primordial magnetic seed fields. 
In this paper, we now specialize to discuss the QCD 
domain walls presented in |po|| showing that they may 
be able to generate sufficiently large magnetic fields to 
seed galactic dynamos and possibly observe in the extra- 
galactic medium. 



C. Strong CP Violations in Domain Walls 

We present a brief summary here of the essence of 
strong CP violation to explain how axion and QCD do- 
main walls may satisfy the criteria discussed above. The 
most general form for the fundamental QCD Lagrangian 
is known to contain the following term related to the 
anomaly: 



32^2 ^M-^ 



(3) 



where 



is the gluon field tensor and G"^" = 
^e^'^ C^p^ is its dual. This term is odd under the discrete 
symmetry CP, thus, if 6 is non-zero, then the strong in- 
teraction should violate CP. Experiments, however, have 
placed tight limits 16*1 < 10"^. The contributions to the 
final 9 arise from several sources, and it remarkable that 
these seem to exactly cancel. The origin of this cancella- 
tion is known as the strong CP problem. 

One solution is to promote 9 from a parameter to a dy- 

. The idea is that, 



namical field called the axion |34 



prior to the QCD phase transition, the axion field is mass- 
less and 9 can take on any value. After the transition, the 



axion acquires a mass and sits in a potential with a mini- 
mum energy where — Q. The axion field thus relaxes to 
the minimum restoring CP conservation today. To date, 
axions have not been detected, however, there is an al- 
lowed region consistent with experimental, astrophysical 
and cosmological constraints: Thus, the so-called invisi- 
ble axion ^ , which is very light: ~ 10~^-10~'^ 
eV, may still resolve the strong CP problem. In addition, 
axions of this mass are very strong cold dark candidates 
(see for example the recent review ||3^). As mentioned, 
axions provided a nice mechanism for generating domain 
walls, but because the axion must be so light, there is no 
way for such structures to efficiently interact with nucle- 
ons. For a good reviews of the strong CP problem and 
the role of axions, see |39, 40, 



In any case, we assume that some method exists to 
solve the strong CP problem. What is important about 
the 9 parameter is that, in the low-energy limit, it only 
appears in the combinations {9 + (p) (Equation ^ and 
{9 + (p — a), where is the dynamical field related to 
the 77' meson and a is the axion field. Thus, even when 
9 — 0, CP will be violated in strong interactions in a 
domain wall background where (/) or — a is non-zero 
over a macroscopically large region. Hence, QCD and 
axion domain walls induce strong CP violations over their 
central regions. 

One of the consequences of this strong CP violation is 
that nucleons have an induced electric dipole moment as 
well as a magnetic dipole moment"* pl| ]. We summarize 
those results here. 



In the chiral limit nia 



for small 9 



47r2 



rriN 



(4) 



where q-kMN is the strong ttNN coupling constant and 
QttNN is the CP odd ttNN coupling constant which was 
estimated to be (jTrNN ~ 0.0416*1. In these formulae the 9 
parameter should be treated as the singlet (p domain wall 
solution (t>{z) with nontrivial z dependence. From these 
formulae one can compute the following relation 



9{z) ~ 0.1. 



(5) 



Thus, for all nucleons, including the neutron, both the 
electric and magnetic dipole moments arc non-zero and 
of the same order in the domain wall background when 
9iz) = <p{z) ^ 1. 



V. QCD DOMAIN WALLS 



We saw in section (IVB) that several types of domain 
walls might act as sources for seed fields. To be concrete. 



"* Normally the electric dipole moment is suppressed to the same 
order as 6 in that CP is conserved. 
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we shall now restrict our attention to QCD domain walls 
to show how domain walls might produce magnetic seed 
fields. In this section, we shall present a short review of 
the results presented in simplifying the model for 
presentation. 

To describe these walls, we consider the low-energy 
effective theory of QCD including the pions and the rj' 
singlet field. The rj' field is not as light but is the source 
of the physics behind the QCD domain walls. The pions 
and rj' enter the Lagrangian through the matrix repre- 
sentation 



U 



exp 



iV2—— + i- 

Jit 



(6) 



where 7r° are the N"^ — 1 pseudo-Goldstone fields, A° are 
the Gell-Mann matrices for SU{Nf) and rj' is the singlet 
field. From now on, we limit ourselves to the simplest 
case of one fiavour Nf — 1 which contains only the rj field 
but captures all the relevant physics. The Nf = 2 case is 
presented in . Although the models are quantitatively 
different, the phenomena described by both is the same. 
In this model, we see that (H) reduces to a single complex 
phase 



U = e"^, <^=^. (7) 

The effective Lagrangian density then reduces to 

>c = ^(a^0)'-n0) (8) 



with the effective potential 
V{(j)) = — min I M cos 4> + E cos ( — ^ 



(9) 



which was first introduced in The minimization on 
the right comes from choosing the lowest energy branch of 
the multi-valued potential. Details about this potential 
are discussed in the original paper |^ but several points 
will be made here. All dimensionful parameters are ex- 
pressed in terms of the QCD chiral and gluon vacuum 
condensates: M = diag(m^|(g*(7')|) is the mass matrix 
and E — {bas/{32TT)G^) is the vacuum energy. These 
are well known numerically: niq ~ 5 MeV (see for exam- 
ple PI), (gg) (240 MeV)3 @, (as/7rG2) ~ 0.012 

GeV^lll and b = llNc/S - 2Njf3 is the first term in 
the beta function. 

This potential correctly reproduces the Di Vecchia- 
Veneziano-Witten effective chiral Lagrangian in the large 
Nc limit 47 1, it reproduces the anomalous conformal 
and chiral Ward identities of QCD, and it reproduces the 
known dependence in for small angles |46, It also 
exhibits the correct 2'k periodicity in 9. This periodic- 
ity is the most important property of the potential and is 
the reason that QCD domain walls form: The qualitative 
results do not depend on the exact form (cosine) of the 



potential. Rather, the domain walls form naturally be- 
cause of the 2tt periodicity {9^6 + 2tt) which represents 
the discrete nature of the ground state symmetries. It is 
exactly these symmetries which leads to the existence of 
axion domain walls when 9 is promoted to a dynamical 
axion field Isil, Is^, li^, Egl, |50tl . 



As described in section ( I V C ) , we see that the singlet 
U{1) field (j) occurs in the same place as the CP violating 
9 parameter. Thus, even though to a high degree of pre- 
cision, = 0, in the macroscopic regions where (77') ^ 
is non-zero, there will be CP violating physics. 

We see that the potential(||) has ground states charac- 
terized by 



4>Q = 27rn, 

where n is an integer, with vacuum energy 
E. Expanding about the minimum = c/jq + 
the mass of the field^ 



(10) 



n = -M- 

we find 



(11) 



The most important point to realize is that all of the 
ground states (T^ in fact represent the same physical 
state U = 1. Thus, it is possible for the cf) field to make 
a transition 27rn — > 27rm for different integers n and to. 
Within this model (^), where all heavy degrees of free- 
dom have been integrated out, these transitions are ab- 
solutely stable and represent the domain walls. When 
one includes the effects of the heavier degrees of free- 
dom, however, we find that the walls are unstable on the 
quantum level. This is described in detail in and 
briefly reviewed in the next section. 



A. Domain Wall Solutions 

To study the structure of the domain wall we look at a 
simplified model where one half of the universe is in one 
ground state and the other half is in another. The fields 
will orient themselves in such a way as to minimize the 
energy density in space, forming a domain wall between 
the two regions. In this model, the domain walls are 
planar and we shall neglect the x and y dimensions: A 
complete description of this wall is given by specifying 
the boundary conditions and by specifying how the fields 
vary along z. 

We present here the two basic domain wall solutions. 
These are characterized by interpolations from the state 
(h = Oto: 



Soliton: (p = 27r, 



In the more general case of N t quarks with equal masses, the 
right-hand side of Equation (hlj) should be multiplied by a factor 



TV, 
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Anti-soliton: 



-27r. 



It is possible to consider transitions between further 
states (i.e. 2TTn) but these can be thought of as 
multiple domain walls. They also have higher energies, 
are less stable, and are thus less important for our discus- 
sion. To gain an understanding of the structure of the 
domain walls we look for the solution which minimizes 
the energy density of the domain wall. The energy den- 
sity (wall tension) per unit area is given by the following 
expression 



(12) 



where the first term is the kinetic contribution to the 
energy and the last term is the potential. Here, a dot 
signifies differentiation with respect to z: <i = . 

To minimize the wall tension ( |l2|) , we can use the stan- 
dard variational principle to arrive at the following equa- 
tions of motion for the domain wall solutions: 



E 



AM 



MNr 



(13) 



Again, the last term of Equation ( |l3| ) should be under- 
stood as the lowest branch of a multi-valued function as 
described by Equation (||). 

The general analytical solution of Equations (|3|) is 
not enlightening and we present the numerical solution in 
Fig. 0. In order to gain an intuitive understanding of this 
wall, we examine the solution in the limit M <C E/N"^ 
(physically, when Nf > 1, this is the limit <C m^y)- 
In this case, the last term of ( p^ dominates. Thus, the 
structure of the (j) field is governed by the differential 
equation: 



AE 



■ sm - 



(14) 



Now, there is the issue of the cusp singularity when (fi — tt 
because we change from one branch of the potential to 
another (see Equation (^.) By definition, we keep the 
lowest energy branch, such that the right-hand side of 
Equation ([l4| ) is understood to be the function sin((/)/7Vc) 
for < < TT and sin([0- 27r]/iVc) for tt < (/) < 27r. How- 
ever, we notice that the equations of motion are symmet- 
ric with respect to the center of the wall (which we take 
as z = 0), hence cj) = tt only at the center of the wall and 
not before, so we can simply look at half of the domain, 
z £ (—00,0], with boundary conditions 4>{—oo) = at 
z = — oo and (j){0) = tt at z = 0. The rest of the solution 
will be symmetric with (f) = 2^ at z = +oo. 

Equation (|lj) with the boundary conditions above has 
the solution 



which is a good approximation of the solution to Equa- 
tion (^ when M < E/N^. Here, the scale of the waff 
is set by the parameter ^: 



(16) 



2VE 

— — , hm ii = mr,' 



which is the inverse width of the wall and which is equal 
to the field mass m,,' in the chiral limit — > (see 
Equation (pi])). Thus, we see that, indeed, the QCD 
domain walls have a QCD scale. 

Solution ( |r^ ) describes the soliton. The anti-soliton 
can be found by taking z ~z: thus, we have the tran- 
sition soliton^anti-soliton under the discrete CP sym- 
metry. The numerical solution for the (j) field is shown in 
Fig. |l|. It turns out that the approximation is reasonable 
even in the physical case where N^M/E 10~^. 

The wall surface tension defined by Equation (|l^) and 
can be easily calculated analytically in the chiral limit 
when the analytical solution is known and is given by 
Equation (15). Simple calculations leads to the following 
result:^ 



a ^ 4NJ, 



(17) 



In the case when ruq ^ 0, a is numerically close to the 
estimate (^. 







FIG. 1: Basic form of the QCD domain wall (soliton). The 
analytic approximation ( [l5| ) is plotted as a dotted line to show 
the good agreement. We have taken A'^c = 3 here. Notice that 
the wall thickness is set by the parameter jj.. 



(^) 



271 - ANr tan" 



e-M^an^^ 



z < 0, 
z > 0. 



(15) 



" In the general case of Nf quarks of equal mass, the right-hand 
side of Equation should be multiplied by the factor 1/ ^ Nf. 



In this case, it reduces to the result cited in 



when TV 
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B. Domain Wall Decay 

Finally, we note that these domain walls are not stable: 
as mentioned earlier, the vacuum states ( p^ ) represent 
the same physical state. When one includes the heavier 
gluonic degrees of freedom, it becomes possible for the 
fields to "unwind" through this extra degree of freedom. 
Classically this is not allowed because the heavy degrees 
of freedom are constrained by a large potential barrier, 
but it is still possible for the field to tunnel through this 
barrier forming a hole in the domain wall. Once a large 
enough hole is formed, it will expand and consume the 
domain wall. This process is called "nucleation" and is 
similar to the mechanism consuming iV = 1 axion domain 
walls 111 H, |8[|l). 

In ]20| , we estimate the lifetime of these domain walls 
borrowing the same methods used to estimate the life- 
time of axion domain walls in the A'^ = 1 axion models 
||9|, H, |8[ l^. We should point out one major differ- 
ence, however, between the iV = 1 axion model and our 
model. In the axion model, prior to domain wall forma- 
tion there is a phase where stable axion strings can form. 
When the domain walls form, these strings bound the 
domain walls and thus the walls start to decay from the 
outset greatly reducing their lifetime (see for a nice 
discussion). In the case of QCD domain walls, strings 
are not stable objects and do not form before the walls. 
Thus, the only way for the walls to decay is through the 
nucleation process we are about to describe. This greatly 
suppresses the decay rate and is the source of the long 
lifetime for the walls. 

The tunneling probability can be estimated by comput- 
ing the action S'o of an instanton solution of the Euclidean 
(imaginary time, t — it) field equations, approaching the 
unperturbed wall solution at r — > ±oo. In this case the 
probability P of creating a hole is proportional to the 
factor 



P ■ 



'So 



(18) 



where S'o is the classical instanton action. 

If the radius Rc of the nucleating hole is much greater 
than the wall thickness, we can use the thin-string and 
thin-wall approximation. In this case, the action for the 
string and for the wall are proportional to the correspond- 
ing world-sheet areas [pl[. 



n 4:71 , 

So = ^T^R^a - -jR CF- 



(19) 



Here a is the wall tension (|17|), and a ^ \/2E is the 
string tension which we estimate based on dimensional 
arguments. The string tension a tries to close the hole 
while the wall tension a tries to widen the hole. Mini- 
mizing (nS) with respect to R we find that 



If a hole forms with radius R > Rc then the hole will ex- 
pands with time as x^+y"^ — R'^+t'^, rapidly approaching 
the speed of light and consuming the domain wall. 
Inserting numerical values for the phenomenological 



relevant case Nf — 2 we find that |20 

a - (0.28 GeV)2, a ^ (200 MeV)^, S'o - 120. (21) 

What is important is that So is numerically large, and 
hence the lifetime is much larger than the QCD scale 
because of the huge tunneling suppression e~^° ~ 10~^^. 
A more complete analysis is presented in [^o| where we 
estimate the lifetime of the walls to be of the order 



10"^ s 



(22) 



even though the walls are governed by the microscopic 
QCD scale. This result should be interpreted with some 
caution: in the low energy regime, we do not have very 
good control over the quantitative physics.^ Arguments 
presented in ||2^ show, however, that it is at least pos- 
sible for domain walls of purely QCD origin to live for 
macroscopically large lifetimes. 

To summarize, we have a QCD domain wall with all of 
the properties required to generate magnetic fields: 

1. The walls form shortly after the QCD phase tran- 
sition and attain Hubble-scale correlations through 
the Kibble mechanism. 

2. The QCD domain walls have a structure on the 
scale of m^} ^ ^qcd ^^^^ they can efficiently 
interact with nucleons and other QCD matter. 

3. The transition in the singlet rj' field produces an 
environment near the wall where the effective 6 pa- 
rameter is non-zero. Thus, across the wall, there is 
maximal CP violation. In such an environment, it 
is known that the electric and magnetic dipole mo- 
ments of the nucleons are of the same order . 

4. The strong CP violation also provides a mechanism 
for generating helicity on a Hubble scale^ by align- 
ing both the electric and magnetic dipole moments 
along the domain wall. 

5. The decay mechanism described in pO| renders the 
QCD domain walls unstable such that the walls 
themselves do not pose a cosmological problem. 
However, the suppression in the decay mechanism 
due to quantum tunneling might extend the life- 
time of the walls to a macroscopic scale ( p^ ) which 
is long enough to generate the required electromag- 
netic turbulence as we shall show. 



Rc — — , 
a 



(20) 



^ It is possible to regain control of the calc ulation in the high 
density limit. See the discussion in section ( VIII A| ). 

* To be precise, the domain walls separate the helicity into Hubble 
size regions to that globally the total helicity is zero, but within 
Hubble scale regions, the helicity is maximal and correlated with 
the same sign. 
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VI. ALIGNMENT OF SPINS IN THE DOMAIN 
WALL 



Now we shall show that the domain walls indeed do 
acquire a magnetization and present a simplified method 
for estimating the magnitudes of bulk properties on the 
domain wall. This method makes the approximation that 
the domain wall is flat and that translational and rota- 
tional symmetries are preserved in the plane of the wall 
(which we take to be the x-y plane.) These approxima- 
tions are valid in the case of domain walls whose cur- 
vature is large in comparison to the length scale of the 
pertinent physics. 

Once this approximation is made, we can reformulate 
the problem in 1 + 1 dimensions {z and t) and calcu- 
late the density of the desired bulk properties along the 
domain wall. To regain the full four-dimensional bulk 
properties, we must estimate the density of the particles 
in the x-y plane to obtain the appropriate density and 
degeneracy factors for the bulk density. Thus, the final 
results are not independent of physics in the x-y plane, 
but rather, these effects are accounted for only through 
the degeneracy factors. 

We proceed to demonstrate this technique by calculat- 
ing the alignment of fermionic spins along the wall. We 
take the standard form for the interaction between the 
pseudo-scalar r]' field and the nucleons which respects all 
relevant symmetries: 



£4 = * (i^- 



(23) 



Here 4> = (j){z) characterizes our domain wall solution as 
expressed in Equation (^5|) and mi^ is the nucleon mass. 
For our approximations, we assume that fluctuations in 
the nucleon field do not affect the domain walls and, 
thus, treat the domain walls as a background field^. The 
strategy is to break (p3|) into two 1-1-1 dimensional com- 



ponents by setting dx = dy = (this is the approxima- 
tion that the physics in the z direction decouples from 
the physics in the x-y plane) and then by manipulating 
the system of equations that result. 

First, we introduce the following chiral components of 



the Dirac spinors 



1 
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Xi 



^--Vs (I) ' 




V2 l*H 



(24) 



(25) 



where S is the area of the wall. This normalization factor 
cancels the degeneracy factor proportional to 5 added in 
the next subsection. Now we re-express ( p3| ) by noting 
that 7I = /: 



do -ajdj 
(Tjdj -do 



ruN 



cos((/)) i sm{(j>) 
i sin(0) cos((/)) 



4-. 



The associated Dirac equation is 



do 



-do 



^ cos(0) isin(0) 
^ ' isin(0) cos(0) 



* = 0. 



This is equivalent to the coupled system: 

2i{do + cr.a,)^'- = 2mAre"^*+, (26a) 
2i{do - (7^90*+ = 2mAre"'"^*_. (26b) 

Now, we decouple the z coordinates from x and y by 
setting dx = dy = 0: 



■ fdt + 9. 







dt-d, 
dt - d, 



a.-Ka.r+ = --^""*- 



(27a) 
(27b) 



Both these equations are diagonal. Thus, we see that 
the top components and bottom components of mix 
independently. 



-mjve"^ iidt + d,)\ (xi 

i{dt-d,) -mNC-^y y^i 

-ruNe'^ iidt-d,)\ fx2\ _ 
i{dt + d,) -ruNe-'^ U2 



= 0, 
0. 



(28a) 
(28b) 



Remember that we are looking for a two-dimensional 
Dirac equation, thus we want the kinetic terms to look 



A full account would take into account the effects of this back- 
reaction. We expect that such back-reactions would affect the 
potential (Bl) by altering the form of the last term and possibly 
adding higher order corrections. This may affect the magnitudes 
of some of the estimates, but would certainly not alter the topol- 
ogy of the fields and thus the domain walls would still form with 
a similar structure. Quantitatively this would alter the results, 
but not the order of magnitude. 



We are using the standard representation here: 

^» = (0 -l) ' = (A ) ■ ^5 = (/ J) ' 

= (? 0) ' = (° 7) ' = (0 -1) • 
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the same. For this reason we should flip the rows and 
columns of the second equation. Doing this and defining 
the two two-dimensional spinors 



(1) 



(2) 



(29) 



the equations have the following structure: 

(i7^a^ - mjve+'^^= = (30a) 
{trOf, - mAre-^"^^=)^(2) = (30b) 

where the index fi £ {t,z}, the Lorentz signature is 
(1,-1) and we define the following two-dimensional ver- 
sion of the gamma matrices: 



7t = CTi 



75 = 0-3. 



These satisfy the proper two-dimensional relationships 
75 = 7t7z and = g^i. + e^,.75- We can re- 

produce equation ( |30| ) from the following effective two- 
dimensional Lagrangian density, 



A =^(1) (ird^ - mAre+^*^^) 



(31) 



where two different species of fermion with opposite chi- 
ral charge interact with the domain wall background 
0(z). Note that, due to the normalization factor l/^/S 
we introduced above, the two-dimensional fields \E'(i) have 
the correct canonical dimension 1/2. 

We have thus successfully reduced our problem to a 
two-dimensional fermionic system. It is known that for 
several systems in 1 -I- 1 dimensions, the fermionic repre- 
sentation is is equivalent to a 1 -I- 1 dimensional bosonic 
system through the following equivalences^^ M, p5[: 



*b)n5*0) 



TT 

'^jbcos(2^/TT9j), 
-^sm{2y/Tr9j). 



(32a) 

(32b) 

(32c) 
(32d) 



After making these replacements, we are left with the 
following two-dimensional bosonic effective Lagrangian 
density describing the two fields 0i and 62 in the domain 
wall background (j){z) 



C 



1 



1 



(33) 



The constant fj, in the last two equations is a scale parameter 
of order mjv • The exact coefficient of this term depends on the 
model an is only known for exactly solvable systems but in all 
cases, is of order unity. This technique is well-known to the 
condensed matter and particles physics communities. See for 
example p3| . 



where the effective potential is 

U{0i,O2) = -niNU [cos(2^/¥6'l 



cos(2\/7r^2 



(34) 

The next approximation that we make is to neglect the 
dynamics of the 0k fields: we assume that they relax 
slowly in the domain wall background such that their dy- 
namics do not contribute appreciably to the final state^^ 
which minimizes the potential (|3^). The classical mini- 
mizing solution is thus 



We are now ready to show that the domain walls align 
the spins of the fermions. The relevant spin operator is 
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*tg^ = *775* = 1'^(T*+ + *^(f*_ (36) 

Let us consider the z component of the spin. We then 
have 

= ^(^(1)7*^(1) -*(2)7t*(2)) (37) 

and so we see that the four-dimensional spin operator 
= 707x75 is expressed in terms of a pair of two- 
dimensional fermion charge operators. We can calculate 
the expectation value of the spin operator in the domain 
wall background using this two-dimensional correspon- 
dence (^^ and the bosonic representation of the fermions 

= |(^'(i)7t*(i) - 'I'(2)7t^'(2)) 
^^d.{02-0i). (38) 

Finally, we use our minimizing bosonic solution ( |35| ) 
to obtain the following four-dimensional average spin 
aligned along the domain wall: 



1 d4>{z) 



(39) 



Sti dz ■ 

We will also need the following matrix elements later on: 

-i(xki-xk2+elxi-ek2) 

= i(*(l)*(l))-i(*(2)*(2)> 

= (40) 



This is the same adiabatic approximation used by Goldstone and 
Wilczek H. 

Here we use the convention that cTij = ^[ji,"/]], thus the spin 
operator Ej, = ^eij^aij. In matrix form with the standard rep- 



resentation, this becomes: 



-[ai,(Tj 








o-fc 

CTfc 



4 \ u — L"n'^J 

In terms of the gamma matrices, this is Ej. = 7o7fc75 
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and 



^1X2 

(i)i75*(i)> 



-i(T3 
i(T3 

xki - xk2) 



(2) 



f sin(0) 



(41) 



Remember that we have restricted ourselves to a 1 + 1 
dimensional theory. We must now estimate the density 
and degeneracy of the nucleons along the wall so we can 
obtain a true 1 + 3 dimensional estimate of the spin den- 
sity. 



A. Fermion Degeneracy in the Domain Wall 

We have assumed that locally the domain walls have 
only a spatial z dependence. There is still a two- 
dimensional translational and rotational symmetry in the 
x~y plane. These translational degrees of freedom im- 
ply that momentum in the plane is conserved and hence 
we can treat the neglected degrees of freedom for the 
fermions as free. The degeneracy in a region of area S 
will simply be a sum over these degrees of freedom with 
a discrete factor g = 4 = 2 x 2 for spin and isospin de- 
generacy 



^ ' ^ 27r 2tt 

\\<Pf 



47r 



(42) 



Here we estimate the Fermi momentum pp « Aqcd by 
the thermal scale of the fermions"'^'' and assume that the 
Fermi sea is filled^^. 

This completes our estimate of the induced spin along 
the domain wall in a small region S. Combining our esti- 
mate of the spin ( |39| ) from the bosonization scheme with 
the fermion degeneracy ( |4^ ) we obtain the spin density 
along the wall: 



A 



4D 



QCD d(f){z) 

dz 



(43) 



VII. GENERATION OF ELECTROMAGNETIC 
FIELD 

Once the spins are aligned, the nucleon electric and 
dipole moments interact with the electromagnetic fields 
F^n through the interaction 



1 



§(iD^)(ii?^)*. 

(44) 



Here the nucleons have both electric and magnetic dipole 
moments c?,p ~ //.p respectively (||). 

Now, we make the approximation again that the nucle- 
ons align independently of the electromagnetic field, and 
we treat the nucleon field \1/ as a background. The situ- 
ation is a field of dipoles aligned along the domain wall. 
The net fields generated by surface of area willed with 
a constant density of aligned dipoles is proportional to 
since the dipoles tend to cancel. For a perfectly flat 
domain wall of infinite extent, ^ — *■ cxd and thus no net 
field would remain as pointed out in The QCD do- 
main walls, however, are far from flat: the walls have 
many wiggles and high frequency dynamics excitations. 
Thus, the fields generated by the dipoles will not cancel 
on the domain wall, but will be suppressed by a factor 
of (^Aqcd)^^ where ^ is an effective correlation length 
that depends on the dynamics of the domain walls. As 
an upper bound, the extent of the domain walls is lim- 
ited by the Hubble scale. Typically, domain walls remain 
space filling, thus we expect AqJ,j-, < ^ ^ Hubble scale. 
Unfortunately, we presently cannot make a tighter bound 
on ^, however, we shall see that, even in the worst case, 
this mechanism can at least generate feasible seed fields 
for galactic dynamos to amplify. 



The result of sections (VI) and (VIA) is a method for 
estimating the strengths of various sources in the domain 
walls. We now need to couple these to the generation of 
electromagnetic turbulence. To do this properly requires 
the solution to Maxwell's equations as coupled to the 
sources in (^). This is difficult, though no doubt im- 
portant for accurate numerical estimates, and so for an 
order of magnitude estimate we consider a dimensional 
estimate considering the sources as a set of dipoles sitting 
in the domain walls. The spacing between the dipoles is 
set by the QCD scale A^L^ , and the strengths of the field 
can be estimated from (|44|) using dimensional arguments: 



As a check, note that, the dimension here is 3 and, the 
result does not depend on the normalization factor S. 



We are treating the fermions along the wall as a massless, two- 
dimensional Fermi-gas. 

Voloshin obtained a similar estimate |23| . Furthermore, he es- 
timates that there are sufficiently many fermions in the Hubble 
volume to diffuse into the domain wall potential justifying our 
assumption about the filled Fermi sea. 



1 



CAqcd 



(d*(*cr^^75*) • (45) 



This includes the dipole suppression discussed above. 
From (|^, (|l]), (H) and (^ we arrive at the follow- 
ing estimates for the average electric and magnetic fields 



The density of the dipoles is governed by the QCD scale Aqqj-,. 
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(correlated on the Hubble scale) which includes the de- 
generacy factors (remember that fj, « rripf): 



\{E.)\-\{B.) 



0.5 eA^cD _ 10'^ G 



(46) 



This method of estimating the electric and magnetic 
fields produced is extremely crude: we have not solved 
Maxwell's equations, we have not taken back reactions 
into account and we have not fully accounted for the 
motion and geometry of the domain walls. Never the 
less, we expect that the estimates (|6|) to be valid as an 
order of magnitude estimate for the field strengths. The 
approximations we have made and effects that we have 
neglected will be discussed in section ( VIII A ) . Thus, we 
have the approximations for the fields (^6|) which, along 
with (ph, justifies the estimate (pi). 



which affects the magnitude of the fields plays no role in 
disturbing this correlation. 

As we mentioned, eventually, the electric field will be 
screened. The timescale for this is set by the plasma fre- 
quency for the electrons (protons will screen much more 
slowly) LOp which turns out to be numerically close to 
Aqcd near the QCD phase transition. The nucleons, 
however, also align on a similar timescale Ag^^, and 
the helicity is generated on this scale too, so the electric 
screening will not qualitatively affect the mechanism. 
Finally, we note that the turbulence requires a seed which 
remains in a local region for a timescale set by the con- 
ductivity 1^ cr - cT/e^ ~ Aqcd where for T = 100 
MeV, c « 0.07 and is smaller for higher T. Thus, even 
if the domain walls move at close to the speed of light 
(due to vibrations), there is still enough time to generate 
turbulence throughout the Hubble volume. 



A. Helicity 



VIII. CONCLUSION 



Finally, we note that the turbulence discussed in sec- 
tion (HI) should be highly helical. This helicity arises 
from the fact that both electric and magnetic fields 
are correlated together along the entire domain wall, 
(E) ^ (A) / r where (A) is the vector potential and r is 
a relevant timescale for the electrical field to be screened 
(we expect r ~ Aqqj-, as we discuss below). The mag- 
netic helicity density is thus: 



A B 



.2 A 



QCD 



(47) 



It can be seen from ( |45D that both the electric field and 
magnetic fields have the same structure in the domain 
wall. This implies that domain walls (solitons) have the 
same sign of helicity everywhere. Under CP the wall 
becomes an "anti-wall" (anti-soliton) and the orientation 
of the magnetic field B changes direction: thus anti- walls 
have opposite helicity. 

Note carefully what happens here: The total helicity 
was zero in the quark-gluon-plasma phase and remains 
zero in the whole universe, but the helicity is separated so 
that in one Hubble volume where one domain wall domi- 
nates, the helicity has the same sign. The reason for this 
is that, as the domain walls coalesce, initial perturbations 
cause either a soliton or an anti-soliton to dominate and 
fill the Hubble volume. In the neighboring space, there 
will be other solitons and anti-solitons so that there is an 
equal number of both, but they are separated and this 
spatial separation prevents them from annihilating. This 
is similar to how a particle and anti-particle may be cre- 
ated and then separated so they do not annihilate. In any 
case, the helicity is a pseudo-scalar and thus maintains 
a constant sign everywhere along the domain wall: thus, 
the entire Hubble volume is filled with helicity of the same 
sign. This is the origin of the Hubble scale correlations 
in the helicity and in B^. The correlation parameter ^ 



A. Summary 

Here is a brief summary of what we have done and the 
approximations that were made: 

1) In section ( IV B| ) we outlined the generic features 
that domain wall models should posses if they are to 
successfully generate primordial seeds by the method de- 
scribed in this paper. We proceeded with the QCD do- 
main wall model |20 to estimate the magnitude of the 
effects^^. 

2) Making the approximation of thin flat domain walls 
we proceeded to calculate the averages of quantities like 
{^Hcr^i,"^') by reducing the interaction ( p3|) to a 1 -|- 1 
dimensional system. In this approximation we discarded 
the momenta in the x-y plane to get (p7|). We capture 
the effects of these momenta by degeneracy factors in 



If the screening were more efficient, then one might worry that 
the electric field would not last long enough to generate the he- 
licity. 

Further support for the existence of QCD domain walls comes 
from calculations presented in |S7[ | where it is shown that these 
walls almost certainly exist in the high density regime of QCD. 
Thus, domain walls seem to be important features at high den- 
sity. Furthermore, if one accepts a conjecture on quark-hadron 
continuity at low ternperature with respect to variations in the 
chemical potential fj, |58[ , then one can make the following argu- 
ment: If for large /i domain walls exist, but for low fi they do 
not, then there should be some sort of phase transition as one 
lowers II. Thus, the continuity conjecture supports the existence 
of quasi-stable QCD domain walls at lower densities, at least 
down to the densities of hyper-nuclear matter. Coupled with the 
fact that gluon and quark condensates do not vary much as one 
moves from the domain of hyper-nuclear matter to the low den- 
sity limit, one suspects that the qualitative picture holds even for 
zero density. Different arguments based on large Nc counting, 
also support the existence of the meta-stablo QCD domain walls 
at zero chemical potential ;j[po|. 
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section (VIA). This approximation is valid only if the 
physics in the z direction is independent of motion in 
the x-y plane. This approximation breaks down when 
thermal (or other) fluctuations are large enough that the 
physics in the z-t directions no longer decouple from the 
physics in the x-y plane. 

In making estimates like (^9|), another approximation 
we have made is to ignore back reactions. We have 
treated the domain wall as a static background: in re- 
ality, the presence of fcrmions in the domain wall would 
affect the structure. What we say is that such back reac- 
tions will not change the overall structure or scale of the 
phenomenon, however, it will definitely alter the quanti- 
tative results. Thus, estimates like (^) should only be 
taken as qualitative approximations to the structure in 
the domain wall. A comprehensive analysis would take 
into account the effects of fermions on the domain walls 
through additional interactions to (^). These interac- 
tions, however, would not alter the U{X) nature of the 
r{ field, and thus the basic domain wall structure would 
still be present. Also, it is unlikely that the back-reaction 
of the electromagnetic fields can substantially affect the 
domain wall structure or the alignment of the spins. In- 
deed, the domain walls and the spin alignment are due to 
QCD interactions on the scale Aqcd- Any back-reaction, 
would be suppressed at least by a factor of a ~ 1/137, 
thus, the quantitative results might be altered, but we ex- 
pect the qualitative behaviour and orders of magnitudes 
to be preserved. 

3) The next step was to estimate the strengths of the 
generated fields by using dimensional arguments and con- 
sidering a collection of dipoles aligned in the domain wall 
background arriving at the estimates (^). The actual 
fields generated will be sensitive to the geometry and 
dynamics of the domain walls: this is something that we 
need to understand much better. We have captured these 
effects in the unknown scale length ^ but there is much 
we could understand about this. To study these effects 
we will need to solve Maxwell's equations (|4|), however, 
in the non-trivial geometry of the domain walls we will 
probably have to simulate this. Of particular importance 
is the question: Can the larger fields be generated by the 
motion of the domain walls? 

The interaction of the fields with the plasma is also 
important, as the electric fields may be screened. Simple 
estimates, however, show that the screening timescale is 
at least as long as the other timescales. 

4) We have estimated the scaling of the magnetic tur- 
bulence (^) and showed that the fields generated by do- 
main walls at the QCD phase transition would be of as- 
trophysical interest. The point we make here is that fields 
are generated and correlated on the same length as the 
domain walls. Thus, the domain walls provide a mecha- 
nism for generating Hubble scale correlations. 

Without a better understanding of the dynamics of 
the domain walls, we cannot better estimate the field 
strengths, but it is possible that the generated fields 
are quite large (nanogauss to microgauss scale) even 



without any amplification. It seems that the fields will 
still require some amplification by galactic dynamos, but 
these primordial seeds provide the large scale correlation 
lengths that have been difficult to achieve through other 
mechanisms. 

We have seen that domain walls at the QCD phase 
transition may provide a nice way to resolve some of the 
problems with generating large scale magnetic fields. It 
is important to note that the mechanism described here 
is seated in well-established physics and makes definite 
predictions. The only free "parameter" is a correlation 
length ^ which affects only the strengths of the fields gen- 
erated. This parameter is not really free, but represents 
our lack of understanding of the formation and dynam- 
ics domain walls. As this understanding is improved, the 
scale of this parameter will be fixed, and the method 
will make a definite prediction about the strengths of the 
fields. 

Thus, this mechanism is testable: it already makes a 
definite prediction about the order of the field correla- 
tion lengths: If fields are observed to have much large 
correlations that 500 kpc, then either a more efficient 
inverse cascade mechanism must be discovered or an- 
other method (possibly inflationary) must be considered 
to generate these correlation lengths. 



B. Speculations and Future Directions 

This mechanism has been a first attempt to describe 
astrophysical consequences of recently conjectured QCD 
domain walls These walls may affect several areas 

of astrophysics, both through the field generation mech- 
anism described here and through other effects. We dis- 
cuss several of these applications below. 

The most promising possibility to test the idea of pri- 
mordial seeds is to measure extra-galactic and extra- 
cluster magnetic fields. This may be possible in the near 
future through measurements of anisotropics in the cos- 
mic microwave background (CMB) radiation spectrum 
p9| , |60| or through non-thermal radiation from Compton 
scattering of relativistic electrons off of the CMB [ |6l| . 
Current CMB observations place upper bounds on 
the strengths of primordial fields, which the fields de- 
scribed in this paper respect. There are other limits 
imposed on the strengths of primordial magnetic fields 
from nucleosynthesis production rates, all of which this 
mechanism respect. For a thorough review, see 

There are many possible consequences of primordial 
magnetic fields discussed in j^] , but we point out a partic- 
ularly interesting possibility that might support domain 
wall generated fields. It seems that certain types of field 
structures might explain the apparent violations of the 
GZK cutoff It is likely that primordial fields from 
domain walls may posses a planar structure that could 
assist in explaining these violations. 

Independent of the magnetic fields, QCD domain walls 
might affect early cosmology in another way. As we dis- 
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cussed in the text, baryons are concentrated on the wall 
causing inhomogeneities in the baryon density. While it 
seems that QCD domain walls will decay before T = 1 
MeV, the inhomogeneity in baryon density must be re- 
distributed by diffusion and other dissipative processes 
and thus inhomogeneities might persist that affect nucle- 
osynthesis. 

Recent measurements of the CMB by BOOMERANG 
^ and MAXIMA |||| lead to a value of baryon density 
ils is larger than the value allowed by the conventional 
model. The agreement can be achieved through inho- 
mogeneous big bang nucleosynthesis (IBBN) if regions of 
baryon inhomogeneity are separated by a distance scale 
of about 70 km (at T = 1 MeV), and if these regions 
have a planar structure with high surface to volume ra- 
tio. Models have been proposed whereby such a pla- 
nar structure can occur [ |67| , however, we suggest that 
QCD domain walls might automatically create this kind 
of structure with the appropriate scale if the baryon dif- 
fusion and other dissipative processes are sufficiently slow 
(or if something extends the lifetime of the QCD domain 
walls beyond the estimate (|2^)). Conversely, nucleosyn- 
thesis provides another constraint to check the validity 
of QCD domain wall properties. Additional work is re- 
quired |68j before any definite statements can be made. 

At this point, we would like to comment on a specu- 
lative application of QCD domain walls that may be of 
significant astrophysical interest. Given that QCD do- 
main walls are very stable at high densities , it is not 
inconceivable that the accumulation of baryons along the 
domain wall by a mechanism similar to that described in 
section (VI) may stabilize the the domain walls so that 
they survive much longer than (p^. If this is the case, 
then QCD domain walls might form the seeds for primor- 
dial baryonic compact objects. Such objects would cer- 
tainly affect nucleosynthesis, the cosmic microwave back- 
ground radiation spectrum and possible structure forma- 
tion. The result, however, would be stable, cold, "invis- 
ible" baryonic matter that would contribute to the dark 
matter. In addition, this matter would have QCD scale 
cross-section recently suggested for dark matter to 
explain several discrepancies with the standard cold dark 
matter picture. Besides that, this matter would offer a 
simple explanation of why observations find ^loark ~ 
to within an order of magnitude. This fact is extremely 



difficult to explain in models where a dark matter can- 
didate is related to a physics independent of baryogen- 
esis. Without further calculations, it is not possible to 
make even qualitative predictions about these specula- 
tions, but, as our understanding improves, such a model 
would be able to make definite predictions. Furthermore, 
the physics of QCD domain walls may lend itself to sig- 
natures that can be tested in relativistic heavy ion colli- 
sions. This might provide some concrete foundations for 
QCD phenomena affecting cosmology and astrophysics. 
Likewise, the accuracy of nuclear abundance measure- 
ments and CMB anisotropy measurements could provide 
serious constraints on the behaviour of QCD phenom- 
ena related to domain walls. Thus, while we cannot yet 
seriously advocate the idea of compact primordial bary- 
onic objects forming from this mechanism, one should at 
least keep them in mind when studying these and similar 
processes. 

We would like to close by emphasizing the relation- 
ship that is developing between astrophysics and particle 
physics at the QCD scale. Particle physics provides con- 
crete models for astrophysical phenomena that, as our 
understanding of fundamental physics increases, has def- 
inite predictive power. Thus, there is the potential to 
ground many astrophysical phenomena in well founded, 
testable physics. In return, astrophysical observations 
provide a means for testing particle physics theories un- 
der conditions not possible on earth. We eagerly await 
the fruitful developments that will result from this recip- 
rocal relationship. 
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